Multiscale theory of valley splitting 
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The coupling between z valleys in the conduction band of a Si quantum well arises from phenomena 
occurring within several atoms from the interface, thus ruling out a theoretical description based 
on pure effective mass theory. However, the complexity and size of a realistic device precludes 
an analytical atomistic description. Here, we develop a fully analytical multiscale theory of valley 
coupling, by combining effective mass and tight binding approaches. The results are of particular 
interest for silicon qubits and quantum devices, but also provide insight for GaAs quantum wells. 
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The two-dimensional electron gas formed at a silicon 
heterointerface underpins the modern electronics indus- 
try. But in spite of its ubiquity, the silicon interface 
exhibits phenomena that are not fully understood, and 
cannot be explained by the conventional effective mass 
theory. In the emerging field of nanoelectronics, quan- 
tum degrees of freedom like spin form the basis for novel 
technologies, such as spintronics [1] and quantum com- 
putation 0, [1] ■ An alternative degree of freedom is asso- 
ciated with the low-lying features in the conduction band 
structure, known as valleys, for Si and other indirect gap 
semiconductors [4] . The degeneracy of the low- lying val- 
ley states is 2-fold for both Si[001]/SiO2 and Si[001]/SiGe 
interfaces, due to mass anisotropy and strain effects, re- 
spectively. For quantum devices, where valley physics 
has been studied most extensively [Hj, the valley states 
may either be "frozen out," in favor of the spin degree 
of freedom , or utilized as qubits 7f\ . In either case, a 
complete understanding of the valley physics is essential. 

One main concern is the origin and the magnitude 
of valley coupling, which lifts the valley degeneracy. 
The two approaches previously applied to this prob- 
lem involve continuum theories like effective mass (EM) 
[1) B EqIj El, [l^ [l3l. and atomistic theories like tight 
binding (TB) [IJ, |l5[. While the former provides in- 
tuition because of its analytical nature, it cannot fully 
account for the fundamentally discrete and atomistic na- 
ture of the valley coupling. On the other hand, atom- 
istic theories give an accurate description of the valley 
splitting from first principles, but they cannot provide 
analytical results, except in the simplest geometries 
In this letter, we bridge the gap between microscopic and 
macroscopic theories through a multiscale technique, and 
we provide theoretical justification for an intuitive ex- 
tended EM theory. We show how atomic scale corrections 
near an interface lead to significant improvements in the 
EM description, even for direct gap semiconductors like 
GaAs. 

It is well known that the conventional EM theory of 
electron confinement in a semiconductor crystal breaks 
down near a sharp confining potential, like a quantum 
well For direct gap materials, the resulting errors are 



typically small and may be treated perturbatively in the 
long- wavelength EM theory [l^. For indirect gap ma- 
terials like Si, the atomic scale physics of the interface, 
which is absent from the EM theory, is also responsi- 
ble for valley splitting, corresponding to an energy scale 
1 meV) that is comparable with other energies of 
interest for quantum devices. Perturbative treatments 
should therefore be undertaken with caution. Here, we 
describe a multiscale approach that efficiently captures 
the atomic scale corrections by means of an effective in- 
terface potential, A. The approach is fully analytical, 
and can be applied to general quantum well geometries. 

For concreteness, we consider the case of a 
SiGe/Si/SiGe symmetric square well. The EM wavefunc- 
tion is written as 181 
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n— ±2 



(1) 



where a„ are the valley composition factors (of no im- 
portance here), k±2 = ±fco are the positions of the val- 
ley minima in the Brillouin zone, Uk^{r) are periodic 
Bloch functions, and F{z) is the long wavelength enve- 
lope. Note that we have only included contributions from 
the two z valleys, as appropriate for [001] strained quan- 
tum wells. The short wavelength physics is contained in 
the fast phase oscillations and the Bloch functions. Note 
that the silicon valley minima are located near the Bril- 
louin zone boundaries, with [l3| = 0.82(27r/a) for a 
Si cubic unit cell of width a = 5.43 A, consisting of four 
atomic planes along [001]. 

We have previously argued that the atomic scale 
physics of valley coupling can be incorporated into the 
long wavelength theory by means of a (5-function poten- 
tial at the quantum well interface Hlj]. Similar arguments 
have also been put forth in Refs. [l3, [3], leading to a 
coupled set of envelope equations of the form [3, la, 13] 
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2m; dz^ 
+ ^A(5(z -Zm)-E 
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F{z) = 0. 
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Here, mi = 0.916mo is the longitudinal effective mass, 
V{z) is the vertical confinement from the conduction 
band offset of Vq, and — zLL/2 are the two quan- 
tum well interfaces. The A term couples the z valleys. 
Since V{z) is constant away from the interface, the so- 
lutions for the ground state envelope function are given 
by 



F{z) 



Acos{qz) {\z\ < L/2) 
Be-P'' (|z| > L/2) 



(3) 



Conventionally, the unknown parameters in Eq. ([3]) are 
determined by matching the envelope function and its 
first derivative on either side of the interface [2lj . How- 
ever, the latter matching condition must be modified in 
the presence of a 5-function interface potential. Integrat- 
ing Eq. ^ over an infinitesimal range about the inter- 
face, we obtain the new matching condition, 

F'_,{L/2)-FL{L/2) 
2mi F{L/2) ' ^ ' 

where (P-) correspond to right-hand (left-hand) 
derivatives. Note that a similar discontinuity in F'{z) 
also occurs at any heterojunction with two different ef- 
fective masses [21| . This effect is unrelated to valley cou- 
pling, and we ignore it below. Indeed, for silicon-rich 
SiGe materials, the mass variations are small and incon- 
sequential for our main results. 

We now construct a multiscale theory for A. To be- 
gin, we note that the conventional EM theory remains 
valid and accurate, except within about one atom dis- 
tance from the quantum well interface. In the vicinity of 
the interface, the EM theory should be replaced by an 
atomistic one. The simplest TB theory that can describe 
valley coupling was derived in Ref. [l4|. The model in- 
volves two bands, with nearest and next-nearest neighbor 
tunneling parameters, ti and t2, respectively. An addi- 
tional onsite parameter describes the confinement poten- 
tial V{z) of the quantum well. Although these parame- 
ters may vary with position, depending on the alloy com- 
position, the variations are small for SiGe and we ignore 
them here. The tight-binding coupling parameters are 
then given by [14] 

i2sin^(fcoa/4) = 2h^/mia^, ti = 4t2 cos(/coa/4). (5) 

Near the top interface, the TB Hamiltonian is given by 

/•. \ 



H 



t2 h Q ti t2 

{) t2 h Q ti t2 

Q t2 h Vo h t2 

^2 ii Vo h t2 



V 



(6) 



(Diagonal elements are highlighted in bold, for clarity.) 
The eigenstates of H correspond to vectors of TB coef- 
ficients (. . . , C_2, C_i, Co, Ci, . . . ). Here, we consider a 
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FIG. 1: (Color online) Diflterences between the effective mass 
(EM) and tight binding (TB) wavefunction solutions are 
shown for a GaAs quantum well. (Only the right half of the 
well is shown.) The dashed line shows the conventional EM 
theory, while the solid line includes multiscale corrections. 
The vertical dotted line marks the quantum well boundary. 
Inset; Multiscale ansatz for envelope functions near the inter- 
face. For Si, 7 atoms are needed. For GaAs only 4 atoms are 
used. Open circles correspond to sites in the barrier region. 



quantum well of size L = a(A''-|-l)/2, containing (2A''+1) 
atoms, and centered at atomic position TV = 0. 

To implement a multiscale theory, we require that the 
TB eigenstates match the EM solutions away from the 
interface. From Eq. ([T]), the two lowest-energy TB wave- 
functions can be expressed as 



Q = (-l)^'\/2cos(jA:oa/4)F„ 
S, - (-l)JV2sin(jfcoa/4)L;-, 



(7) 
(8) 



where the actual ground state depends on the width of 
the quantum well Here, {—\)^l\f2 corresponds to 

the Bloch function for the lowest band of the two-band 
model [13|, while the cosine and sine functions describe 
the exponential phase factors in Eq. ([T]), corresponding 
to the cases (a+zjO-^) — (1,±1)/V2, respectively. The 
envelope coefficient Fj — F{ja/A) is determined from 
the correspondence with Eq. ([3]). The physics of valley 
coupling is captured by the sudden change of slope (i.e., 
the kink) in the envelope at the interface. To facilitate 
calculations, we consider the model parameters shown in 
the inset of Fig. 1. Immediately adjacent to the interface, 
the envelope function exhibits a change of slope, with 
F'_ = (4/a)/3 and F[_ = (4/a)7 on either side of the kink, 
and amplitude -Fjv+i = / right at the interface. Thus, 
Fn ^ f + f}, Fn+2 = f -1, and so on. 

Before solving the multiscale theory, it is illuminat- 
ing to gauge the accuracy of the EM ansatz of Eqs. ([7])- 
dH]). This is accomplished in Fig. 2, using the two lowest- 
energy TB eigenstates to infer the envelope function from 
the relation 2F^ = Cj + Sj. The result exhibits resid- 
ual short wavelength structure, arising from the fact that 
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FIG. 2: (Color online) Dimensionless envelope for the ground 
state wavefunction in a 10 nm Si/SiGe quantum well. (Only 
the right half of the quantum well is shown.) The vertical 
dashed line marks the quantum well boundary. The discrete 
points are obtained from TB theory, as described in the text. 
Spurious short wavelength structure is cause by limitations in 
the EM theory. Inside (outside) the quantum well, the blue 
(red) solid lines are fits to Eq. ([3]). Inset: the full, dimension- 
less TB wavefunction, including the fast oscillations. 



the exact k values of the fast oscillations of the two low- 
est eigenstates are nearly (but not quite) identical flE^ . 
This is a signature of the incomplete separation of the 
long and short wavelength physics, and it places a fun- 
damental limit on our ability to match the TB and EM 
theories. In the present work, the spurious "jitter" ev- 
ident in Fig. 2 leads to errors in the evaluation of the 
kink. 

The multiscale theory involves just two equations from 
the full TB Hamiltonian, and we can choose which equa- 
tions to use. The spurious jitter in Fig. 2, could be miti- 
gated by an averaging procedure. Indeed, in this way, we 
obtain excellent agreement with previous estimates of the 
valley splitting [13|. However, such techniques detract 
from the simplicity of the multiscale approach. Here, we 
take a different tack, noting that the alternating behavior 
of the TB coefficients in Fig. 2 can be partially mitigated 
simply by using alternating TB equations. We consider 
the following equations centered symmetrically around 
the interface: 

t2CN-2 + tlCN-l + tiCN+l + t2CN+2 = EtbCn, 

+ tiCN+i + VoCjv+2 (9) 

+tiCN+3 + t2CN+4 = EtbCn+2- 

Note that either the cosine ([7]) or sine ^ functions can 
be used here. 

It is necessary to take into account the curvature of 
the envelope functions in system ([9]), in order to avoid 
unphysical solutions. Near the interface, the cosine en- 
velope in Eq. ([3]) has a vanishing second derivative, but 
the exponential envelope does not. To leading order, we 
can express the latter in terms of parameters / and 7 



as follows: Fn+2 = f ~1, Fn+s = f - "^1 + 7^//, and 
F/v_|_4 = / — 37 + 87^//. Evaluating system we now 
obtain 

f;(L/2)-F:(L/2) ^ (/3-7) 4 
FiL/2) f a 



(10) 
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where we have dropped higher order terms in the small 
parameter a/L. The anticipated linear dependence of A 
on Vq [i3| emerges from Eq. ([4]): 



A ~ Voa/A = (1.36 x 10"^°)Vb, 



(11) 



where A is in units of eVm when Vq is in eV. 

Eq. (jlip is our main result, obtained through a mul- 
tiscale analysis of the kink of the envelope function. We 
can compare this with the apparent kink obtained by fit- 
ting Eq. ^ to the full TB envelope function, as shown 
in Fig. 2. In spite of the spurious jitter, the two esti- 
mates agree to within 20%, in the wide quantum well 
limit. We can also compare Eq. (fTTI) to the estimate 
A ~ (7.2 X 10-")yo obtained in Ref. [ii], by fitting the- 
oretical EM predictions to TB numerical solutions for the 
valley splitting. The latter differs from Eq. (fTTj) by about 
a factor of two, which we attribute to the fundamental 
limitations of the EM ansatz of Eqs. ©-IH]). Neverthe- 
less, it is clear that the EM theory and the multiscale 
analysis, described here, capture the essential physics of 
valley splitting, and enable semi-quantitative predictions. 
Thus justified, the more accurate numerical estimate for 
A in Ref. ;13| forms the basis for a fully quantitative EM 
theory. Further improvements in the EM ansatz and the 
kink analysis should lead to better correspondence be- 
tween the estimates for A. 

We now turn to direct gap materials, such as GaAs. 
Although a sharp confinement potential does not cause 
valley coupling in this case, there are still atomic scale 
corrections to the EM theory. As in the Si case, the cor- 
rections tend to be more significant for narrow quantum 
wells 2M- Foreman has also noted that the corrections 
can be treated within the EM theory by introducing a 
(5-function at the interface 17[. Here, we apply the mul- 



tiscale theory developed for Si to the GaAs quantum well, 
obtaining an analytical expression for the GaAs interface 
potential, A. We also compare the improved wavefunc- 
tion solutions to those obtained from TB theory. 

The simplest TB theory for a single (r)-valley mate- 
rial involves just the nearest-neighbor tunneling param- 
eter ti — ~8h? / m* , where a = 5.64 A is the width of 
the GaAs cubic unit cell. For Ala;Gai_a;As, used in the 
barriers, the effective mass m* depends on the composi- 
tion a; to a greater degree than silicon alloys. Hence, ti 
depends on the atomic position. The onsite parameter 
e{z) — Wfi^ /m*a^ + V(z) also depends on composition. 
However, for the sake of transparency, we will ignore ef- 
fective mass variations, taking e(z) = V{z), and setting 
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m* to an appropriate average of the effective masses near 
the interface. Indeed, m* eventually drops out of the 
leading order expression for A, and a more careful treat- 
ment provides only small corrections. 

The matching condition, Eq. (ji]), also hold for GaAs. 
However, because there is only one valley, and the TB 
model has only one band, the TB envelope function and 
wavefunction are now identical: Cj = Fj . Near the inter- 
face, we parametrize the envelope as shown in the inset 
of Fig. 1. The TB eigenstates are smooth (in contrast 
with Fig. 2), so we may now use adjacent TB equations: 

tiCN-i + IiCn+i = EtbCn, 

IiCn + VoCat+i + tiCN+2 — EtbCn+i- 

In this case, there is no need to consider curvature of the 
wavefunction. Solving for /3 and 7 directly, and noting 
the relation Etb = 2ii + £; between the TB and EM 
energies, we obtain the GaAs result, 

A = a{V^ - E)/A. (13) 

In the physically relevant limit of _E ^ Vb, the GaAs and 
Si interface potentials have identical forms. For GaAs 
materials parameters, we obtain A ~ (1.41 x 1Q~^^)Vq. 

The effect of the interface potential on the wavefunc- 
tion is shown in Fig. 1, where we plot the differences 
between the approximate (multiscale) and exact (TB) 
results for a narrow quantum well. We also show results 
for the conventional (A = 0) EM theory. Deviations from 
the EM theory are small in both cases. However because 
there is no jitter in the GaAs TB envelope function, we 
find that Eq. (fT3|) captures the atomic scale corrections 
with great accuracy. 

In conclusion, we have demonstrated that leading cor- 
rections to the effective mass theory at a sharp quantum 
well boundary arise from the atomic scale physics near 
the interface. The corrections appear as a small kink in 
the envelope function. A multiscale approach, combin- 
ing effective mass and tight binding theories, leads to an 
analytical expression for the effective interface potential 
in silicon, which determines the valley splitting. Simi- 
lar corrections apply to GaAs quantum wells, although 
there is no valley coupling. For other device geometries, 
including graded interfaces and electric fields, the present 
approach remains robust. These situations may also be 
treated by a multiscale analysis. 
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